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Time evolution of entanglement of N quantum dots is analyzed within the spin- 1/2 van der Waals (or 
Lipkin-Meshkov-Glick) XY model. It is shown that, for a single dot initially excited and disentangled 
from the remaining unexcited dots, the maximum bipartite entanglement can be obtained in the 
systems of TV = 2, ■ ■ ■ , 6 dots only. 



Introduction. Entanglement is one of the most pro- 
found concepts of quantum mechanics. But only recently 
it has been recognized as a possible resource for quantum 
information processing including quantum computation, 
quantum teleportation, quantum dense coding, or certain 
types of quantum key distributions jjj . 

To describe a bipartite entanglement of quantum-dot 
systems in a pure state, given by pab = {\^P)(^\)ab, we 
apply the von Neumann reduced entropy defined as 0] 

E[pab] = -Tr{/5A log 2 Pa} = -Tr{p B log 2 p B } 

where pa = Ttb{pab} and p B = Tr A {Pab} are the 
duced density matrices. 
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Model. We analyze entanglement in a system of N 
two-level quantum dots described by the Pauli spin op- 
erators (a*, 5*, a*) for each (n = 1, • • • , JV) lattice site 
by assuming that each dot is coupled to all others with 
the same strength independent of separation distance and 
nature. Such a model is referred to as the spin van der 
Waals (SVW) @, Lipkin-Meshkov-Glick g, equivalent- 
neighbor or infinitely-coordinated spin system. Some as- 
pects of entanglement in the SVW systems in different 
contexts have already been analyzed in Refs. (5[||. 

In particular, the quantum SVW system in the XY-likc 
regime can be defined by the interaction Hamiltonian || : 



The creation and annihilation spin operators are given 



by ol 



&n± a%, respectively. Let's assume that the 



initial state describing a system of M (M — 0, • • ■ , N) 
dots excited and N — M dots in the ground state is given 
as |^(0)) = {\1)® m }a{\Q)® {n - m) }b- Then, the solution 
of the Schrodinger equation of motion for the SVW model 
reads as 
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where M' = min(M, N — M). The states in curly brack- 
ets, {|!)®(^-™)|0)® m }, denote a sum of all possible 



[i.e., ( m )] M-dot states with the excitation number 
(M — rn). The time-dependent superposition coefficients, 
C% M (t), are given by 
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C™(t) = -£b™exp{i[n(N- 
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where are binomial coefficients. In the following we 
will focus on the evolution of the system of N dots in 
which a single dot (M — 1) is initially excited. 
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FIG. 1. Evolution of the entanglement, E NM (solid 
curves), and the Schmidt coefficients: pNM (dot-dashed 
curves) and P^ M (dashed curves) for systems of N = 2, • • • , 8 
quantum dots with only one (M = 1) of them initially excited. 
It is seen that the maximum entanglement corresponds to the 
Schmidt coefficients mutually equal or, in general, the least 
different. 
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(a) N=2..6 M=1 (b) N=7..12 M=1 
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FIG. 2. Entanglement as a function of the rescaled time 
Nhd. It is seen that one ebit of entanglement can be achieved 
in systems with the total number of dots ranging from two 
(thickest curve in fig. (a)) up to six (thinnest curve in fig. 
(a)) only. For AT > 6, the maximum entanglement decreases 
with increasing N . The evolution of the systems with higher 
N is marked by thinner curves. 
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FIG. 3. The maximum values of the entanglement E Nyl 
as a function of the total number N of quantum dots in 
systems with single excitation. On the scale of the figure, 
maxt _E 7,1 (i) = 0.9997 does not differ from unity. However, 
for higher values of Y, a monotonous decrease of the maxi- 
mum entanglement is clearly visible. 



and nt" — jr. We find that the maximum of entangle- 
ment, equal to E N,1 (t') = 1 ebit, can be achieved for the 
evolution time t' for N < 6 only. For N > 6, real solution 
for t' does not exist. This conclusion can also be drawn 
by analyzing Figs. 1 and 2. 

Entanglement for Y > 6 reaches its maximum at the 
evolution times t" . This maximum value, given by 

E N ' l (t") = j%{N 2 log 2 N — (N — 2f \o g2 {N - 2) 

-2(jV-l)log 2 [4(iV-l)]} < 1, 

monotonically decreases with increasing TV as clearly de- 
picted in Figs. 2 and 3. 

Conclusion. We have studied evolution of quan- 
tum dots in the spin van der Waals XY model. We 
have applied its Schmidt decomposition to investigate 
information-theoretic entropic properties and entangle- 
ment of the quantum dot ensembles. Closer analysis of 
our solution shows that the quantum dots evolve peri- 
odically into the bipartite entangled states. We have 
found for the systems of N dots with a single dot initially 
excited and disentangled from the remaining unexcited 
(N — 1) dots that the maximum bipartite entanglement 
can be observed in the systems of N = 2, • • • , 6 dots only. 
It is worth noting that the same "magic" number N = 6 
has been found in a different context of the equivalent- 
neighbor entangled webs in Ref. Q. 
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Entanglement. The Schmidt coefficients for our so- 
lutions are given by H 



They enable the direct calculation of the bipartite entan- 
glement via the Shannon entropy 



E™(t) = E[W))(m\ 
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In the simplest nontrivial case, for M = 1, the Schmidt 
coefficients reduce to 



Pi NA (t)= 4^W(Yi/2) 

and P 2 JV,1 (t) = l-Pjf'^i). The evolution of the Schmidt 
coefficients, P 7 ^ r,1 (t), and the entanglement, E N - 1 (t), are 
depicted in Fig. 1. 

The quantum-dot systems evolve into the maximum 
entangled state for times given by 



= E NM = 2^ sm(NKt) log 2 
with the roots 

Kt' : 



4(JV-1) 



T yCSC 2 (f Kt) - 1 



arccsc(^V 2 ( Ar - !)) 
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